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Quarkonium production and NRQCD matrix elements 
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Most recent calculations of quarkonium production are based on the NRQCD factorization formalism. This 
formalism is reviewed. To make predictions about specific cross section, universal NRQCD matrix elements need 
to be extracted from experiments. Extractions from different experimental situations are compared, with some 
emphasis on the extraction from LEP. 



1. INTRODUCTION 

Quarkonium is an interesting system since the 
mass of the heavy quark is much larger than the 
QCD scale, 



mQ > Aqcd- 



(1) 



One would hope, therefore, that it would be pos- 
sible to calculate reliably in perturbation theory. 
If that were the end of the story, things would be 
simple. The only difficulty, which makes quarko- 
nium very interesting, is that there is another 
natural small scale in the problem, the relative 
velocity of the heavy quarks, v. For charmonium 
v\ ~ 1/4, while for bottomonium v\ ~ 1/10. 

Since we have a multi-scale problem, with the 
hierarchy 



m,Q ^> rriQV ^> rriQV 



A 



QCD, 



it is useful to use an effective field theory. The 
effective field theory appropriate for a system 
with two heavy quarks is Nonrelativistic QCD 
(NRQCD) HI). As the name implies, there is 
a relativistic expansion in v, as well as the usual 
perturbative expansion in a s . 

Before the advent of NRQCD, quarkonium 
production and decay was calculated using the 
so-called "Color-Singlet Model" (CSM). In this 
model, the QQ pair is treated in a color-singlet 
state, in the v — > limit. It was discovered, how- 
ever, that the CSM did not describe quarkonium 
production very well in certain circumstances. 
For example, the CSM prediction for ip' produc- 
tion at the Tevatron is off by about a factor of 



30! As we will see, NRQCD can easily explain 
the discrepancy in this situation ||. 

Whether NRQCD is the correct effective field 
theory to be used for quarkonium is an interest- 
ing and important question. There are some po- 
tential problems comparing NRQCD predictions 
to data, in particular the polarization. NRQCD 
make a definite prediction that at large transverse 
momenta, ip' should be nearly 100% polarized jj). 
However, it appears that the data is not follow- 
ing that trend ||. So it may be that the a s and 
v expansions do not converge well for charmo- 
nium. This is not to say that NRQCD is incor- 
rect. NRQCD is a valid effective field theory, in 
that as the mass toq goes to infinity, the theory 
correctly reproduces full QCD. It may be that v c 
is too large, or tuqv^ too small, for NRQCD to 
be useful in describing charmonium production. 
One could still hope that it will work for bot- 
tomonium. 

How is it possible to test NRQCD? There is the 
polarization discussed above . Another way is 
to compare NRQCD matrix elements extracted 
from different experiments. 

2. FACTORIZATION FORMALISM 

The physical picture of production in NRQCD 
begins with a hard scattering, in which a QQ 
pair are produced with any spin, angular and 
color quantum numbers. This process can be 
calculated in perturbation theory. Then the QQ 
evolve into the final state quarkonium H in some 
non-perturbative fashion. This is encoded in 
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the NRQCD matrix elements, which must be ex- 
tracted from experiment. 

The NRQCD factorization formula is the math- 
ematical realization of the above picture. A gen- 
eral production process for i + j — > H + X can be 
written as 

dcr = ^2 da(ij -> QQ[n] + X)(0 H (n)), (3) 

n 

where da(ij — > QQ[n] + X) is the cross section 
for producing QQ in state n by scattering i and 
j, calculable in a perturbative expansion in a s 
and perhaps convoluted with parton distribution 
functions (PDFs). The long distance matrix el- 
ements encode the hadronization of the heavy 
quarks in state n into the final quarkonium state 
H. The matrix elements can be written as 

ip H {n)) = (0\^T n x \Y,H + X) 

x 

x(£i/ + AVr"V|o) (4) 

x 

The r™ contains Dirac matrices, color matrices 
and derivatives. The NRQCD matrix elements 
scale with a definite power of v, determined by 
T", allowing the truncation of the v expansion. 
Again, we expect the velocity expansion will work 
better in the bottom sector than the charm sec- 
tor. The factorization described above still needs 
to be put on firmer ground. Wc should keep in 
mind that when we use the NRQCD factorization 
formalism to predict production rates, we are not 
only testing the a s and v expansions, but also the 
factorization formalism. 

The scaling is determined by looking at what 
perturbations are necessary to give a non- 
vanishing result for the time ordered product. To 
give non-zero overlap, multipole moment inter- 
actions may need to be inserted. For example, 
the matrix element (0%^ ( 3 Si)) scales as w 4 com- 
pared to the (Of^( 3 Si)), since we need two El 
insertions in the amplitude. The first insertion 
changes L and neutralizes the color. The second 
changes L back to an S wave. Each insertion cost 
a factor of v in the amplitude, or a total of v 4 in 
the rate. The scalings of the most important ma- 
trix elements for ip production are 

(of( 3 so> ~ A (5) 



(OtfSx)) ~ v\ (6) 
{OfCPj)) ~ v\ (7) 
(OfCSo)) ~ v\ (8) 

We can relate the color-singlet matrix elements 
to the wave function at the origin, or derivatives 
of the wave function at the origin, using the vac- 
uum saturation approximation 

(o H (n)) « (o\^r n x \H)(n\x j r n m, (9) 

dropping the sum over X. This allows us to use 
the lattice, data, or models to obtain the color- 
singlet matrix elements. It also means that the 
color-singlet model has been incorporated into 
the NRQCD factorization formalism. 

A quick aside. Many people say the NRQCD 
model, or the "color-octet" model. However, 
NRQCD is not a model. The CSM, on the other 
hand, is a model, since there is no limit in which 
it will reproduce full QCD. 

3. QUARKONIUM AT THE TEVATRON 

As previously mentioned, the CSM fails to de- 
scribe the data at the Tevatron. We will now dis- 
cuss how NRQCD improves the situation. We will 
concentrate on J/ip production. A similar anal- 
ysis can be done for T production. See Refs. 
for recent analyses. 

At large transverse momentum, the leading or- 
der NRQCD prediction in a s and v is O(a^v ). 
This is the same as the leading order CSM predic- 
tion, which grossly underestimates the rate. Even 
the shape of spectrum is wrong, with the theory 
dropping as \/p\ compared to the data which 
falls as l/p\- 

We can go to higher order in the perturbative 
expansion: O(a^.v ). In this case, fragmentation 
processes occur, which, due to some propagators 
being close to on-shell, can have large contribu- 
tions. These fragmentation diagrams dominate 
over the lower order graphs at large p± , with the 
theory now falling as l/pj_, similar to the data. 
However, the normalization of the theory is still 
far below the experiment. 

This is as far as the CSM can go. To improve 
the prediction, NRQCD is needed |y. Instead 
of just including higher order in a s , we can also 
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Table 1 

Table for Mi in units of 10~ 2 GeV 3 extracted from the Tevatron data. The second set of errors, when 

present, are due to scale variation. The first set are statistical, k varies between 3 and 3.5. 

Rel MRSDO MRS(R2) CTEQ2L CTEQ4L GRVLO GRVHO 



6.6 ± 1.5 



10.09 ± 2.07+i'26 " 4.38± 1.15to'?4 3.90 ±1.14^;^ 

10| 1.32 ±0.21 - " 1.44 ±0.21 - - 0.60 ±0.21 



11 1 1.32 ±0.21 - 1.32 ±0.21 - - 0.45 ± 0.09 



12 1 - - 6.52 ±0.67 



include higher order in v: 0(a 3 v 4 ). Now there 
are color-octet fragmentation processes propor- 
tional to the matrix element (Og ( 3 Si)), which 
also scales as l/pj_- Of course, the value of the 
matrix element is unknown, so the fact that we 
can get a good fit really means that the extracted 
matrix element is not abnormally large or small. 

The fragmentation contribution can be approx- 
imated as 

da(pp -^H)=J2 I dz d<J * D i^H, (10) 

i J 

where dii is the cross section to produce an on- 
shell parton i, is the fragmentation func- 

tion for parton i to produce quarkonium state H 
with momentum fraction z, and the sum is over 
partons. The advantage of writing it in this way 
is that it is easy to sum large logs of \ogp±/rriQ 
using Altarelli-Parisi (AP) evolution. 

At lower pj_, fragmentation no longer dom- 
inates, and it is necessary to include all non- 
fragmentation diagrams as well. Contributions 
due to other octet matrix elements, (Og ^Sq)) 
and (Og ( 3 Pj)) are now important. These chan- 
nels both fall as l/p±, making it impossible to 
determine them independently. Instead, a linear 
combination of the two 

Mjf = (OtCS )) + \{OtfP Q )) (11) 

is usually extracted. 

At small p±, higher twist contributions become 
important. These contributions scale as powers 

of A/ ' Jp±_ + mQV n , where the value of n is not 
known. Also at small p±, the rates tend to di- 
verge, due to soft gluon effects. To handle these 
divergences, the low p± region can just be ignored 



or modeled by including initial parton shower- 
ing/intrinsic 

In Tabic [lj different extractions of Mu^ are 
collected. The second set of errors, when present, 
are due to scale variation. The other errors are 
statistical. As can be seen, there is large uncer- 
tainty due to the choice of PDF. At this point, the 
best we can hope to say is that this combination 
of matrix elements is on the order of 

~ few 10~ 2 GeV 3 . (12) 

Table || contains different extractions of 
(Og^( 3 5i)). Again the second set of errors, 
when present, are due to scale variation, while 
the other errors are statistical. Here the largest 
error is not due to the PDFs, but due to the scale 
variation. This is an indication that higher or- 
der perturbative corrections may be very large. 
Again we only know the matrix element is on the 
order of 

(O^CSi)) ~ few 10- 3 GeV 3 . (13) 

4. LEP EXTRACTION 

We would like to do better than the order 
of magnitude extractions discussed above at the 
Tevatron. Since the dominant errors are due to 
scale or PDF variation, it seems difficult make im- 
provements at the Tevatron. We would also like 
to test the formalism, by comparing the matrix 
elements extracted from different experiments. 

The natural place to look for cleaner extrac- 
tions would be at an e + e~ machine, since there is 
no PDF nor initial state gluon radiation to worry 
about. The prediction at CLEO does not have 
a large dependence on the color-octet matrix el- 
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Table 2 

Table for (Og ( 3 S'i)) in units of 10~ 3 GeV 3 extracted from the Tevatron data. The second set of errors, 
when present, are due to scale variation. The first set are statistical. 



Rcf. MRSDO MRS(R2) CTEQ2L CTEQ4L 



GRVLO 



GRVHO 



[8] 

ii 
m 
ii] 

12] 



6.6 ±2.1 

14.0 ± 2.2l J. 3 g 5 
2.1 ±0.5 - 3.3±0.5 

6.8 ±1.6 - 9.6 ±1.5 



10.6 ± 1.4: 



3.94 ±0.63 



11.2 ± 1.4 



-5.6 



3.4 ±0.4 
9.2 ±1.1 



ements away from the endpoint. Therefore LEP 
is the natural choice. 

Formally there are two leading order contri- 
butions in the a s and v expansion, both in the 
singlet channel, of order O(a 2 v ). The contribu- 
tion from gluon radiation in the singlet channel 
Z — > i>gg is suppressed by powers of M 2 /E^ Jl3| . 
The color-singlet charm quark fragmentation pro- 
cess Z — ► ipcc [Q, which has no power suppres- 
sion, dominates over non-fragmentation processes 
for large E^. Light quark octet fragmentation 
(in which the mother parton does not combine 
to form part of the bound state) is naively of or- 
der «gf 4 , down by ~ 1/10 compared to charm 
fragmentation. However, this channel is enhanced 
due to the presence of large logs and a numerical 
factor of five due to the number of possible quarks 
that initiate the process |l5| . The same logs that 
enhance the octet channel also put the conver- 
gence of the perturbative expansion into question. 

The tree-level calculation of the differential 
cross section in the color octet production channel 
p| scales as z — > 1 as 

§^!2?»<or <•*», (») 

leading to large double logs in the total rate. 
Since a s log(M§/M|) « 1.5, we should treat 
a s log(A/|/M|) as order one and resum all pow- 
ers of the large logarithm. With this counting, 
the octet channel is 0(a°u 4 ), on par with the 
singlet fragmentation contribution. More practi- 
cally, the tree-level calculation has a factor of two 
uncertainty associated with the scale at which a s 
is evaluated, since a s {M^)/a s (Mz) ~ 2 (this is 
just a restatement that there is a large logarithm). 
The resummation of the leading logarithms re- 



duces this uncertainty, so the resummation pro- 
cedure is essential from both a practical and a 
formal standpoint. 

To resum the logs, we write the rate in the 
fragmentation limit as 

f z = J^[2C q (^,y)D q ^, Z /y) 

+ C g (n 2 ,y)D g ( f i 2 ,z/y)]. (15) 

We can now use AP to sum the logs of M^/M 2 . 

However, summing the above mentioned logs 
will only yield the correct leading order differen- 
tial rate if z is sufficiently large When z is para- 
metrically small, terms of the form a s log(z)/z 
become just as important. Furthermore, these 
logs will also contribute double logs to the total 
rate given that the lower limit on z is 2M^/Mz- 
This second type of log, due to soft gluon emis- 
sion, is resummed using a formalism familiar from 
discussions of jet multiplicities fig] , where the 
color-coherence of the soft gluon emission is very 
important. 

It is possible to do both of these resummations 
p7|] , see Fig. 1. The rate depends on the linear 
combination 

J2(°s im) ( 3 Si)) x Br[H(m) - J/^X] 

m 

= (1.9 ± 0.5 stat ± l.Otheory) x 10~ 2 GeV 3 (16) 

since the data includes feeddown from higher 
charmonium states. It turns out that this leads to 

1 At very large z ~ 1 we do not get the correct result since 
we have neglected the shape function which enters near 
the endpoint. This shape function has been argued to be 
the solution of the so called "Hera anomaly". It is not 
important here since the rate is very small in this region. 
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an extraction with much smaller theoretical un- 
certainty than the analogous value from the Teva- 
tron 1.4 ± 0.2 stat ± \A t heory For the Tevatron 
value, the extraction from Ref. Q was used for 
the central value. 



singlet decay rate predicts the ratio of the first 
moment over the zeroth moment to be 0.62. Re- 
summation softens the color singlet decay rate, 
but the ratio is still too large, 0.47. The ratio is 
independent of the color singlet matrix element. 
A rigorous extraction of the first moment by the 
experimental groups could provide an extremely 
clean, quantitative test of the NRQCD approach. 
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Figure 1. Differential rate dT/dz as a function of 
z = 2E^ jMz vs data. The dashed line is the sum 
of the tree-level octet and singlet results and the 
solid line is the octet plus the singlet resumma- 
tion. From Ref. lOl 



It is difficult to tell from Fig. 1 if there is much 
improvement in the prediction from this resum- 
mation. To make the comparison more quanti- 
tative, it is useful to take the ratio of the first 
moment to the zeroth moment 



1 



dz z — (Z 

dz 



JH>x) 



(17) 



T(Z -> J/ipX) 

which for the resummed rate is 0.30. The tree 
level rate gives ~ 0.5. A very rough estimate of 
this quantity obtained from the data Q suggests 
a value of 0.26 ± 0.10. This is in sharp contrast to 
the color singlet prediction. The tree-level color 



5. OTHER EXTRACTIONS 

5.1. HERA 

The NRQCD prediction does not seem to fit 
the data at HERA very well. As z — Ej/^/E 7 
goes to the endpoint, z -» 1, the NRQCD pre- 
diction begins to rise. This has been dubbed the 
"HERA anomaly", and has been used to argue 
that NRQCD is incorrect. 

However, it is now known that the NRQCD 
prediction breaks down near the boundary of 
phase space, precisely where there is a problem 
at HERA In this region the velocity 

expansion is breaking down, and a shape func- 
tion must be introduced. In fact, one should 
not compare the prediction with the data above 
z ~ 1 — v 2 « 0.7. If one looks only below this 
point, there is no problem with the theory. Un- 
fortunately, it is not really possible to extract the 
matrix elements from the data in this region. 

5.2. B -> J/ip DECAYS 

For B — > J/ip decays, the first worry is that 
the corrections to factorization will be large since 
there is not a lot of energy in the decay products. 
There are also a s corrections |2l[] and corrections 
in the Heavy Quark Effective Theory l/m& ex- 
pansion J22| , the most important being the Fermi 
motion of the b-quark, which requires the use of 
another shape function. 

Nevertheless, it is possible to obtain the lin- 
ear combination from B decays, with qual- 
itative agreement with the Tevatron extractions 
I2I.I2I 



= 1.5 x 10~ 2 GeV 3 , 
= 2.4 x 10~ 2 GeV 3 . 



(18) 
(19) 
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6. CONCLUSIONS 



E. Braatcn, B. A. Kniehl and J. Lee, hep- 



— NRQCD has been nsext for a nnnrbeT of years 

to make predictions about quarkonium produc- 
tion. At this time, there are consistent extrac- 
tions of the NRQCD matrix elements from a num- 
ber of different experiments. The most analysis 
has been done on hadronic collisions, in particular 
at the Tevatron due to the high statistics. 

LEP is a very clean place to extract the ma- 
trix elements. It is possible obtain smaller theo- 
retical uncertainty compared to the Tevatron ex- 
tractions. Unfortunately the statistics are not as 
good. Nevertheless, due to the large theoretical 
errors in the Tevatron extractions, the LEP ex- 
traction is the best extraction to date. 

It is also possible to test the NRQCD formal- 
ism at LEP by looking at the moments of the 
spectrum. Without the color octet channels in- 
cluded in the NRQCD formalism, the theory and 
experiment would not agree for this test. This 
is one place where NRQCD passes a test. It will 
be interesting to see whether it will be possible 
to resolve the polarization problem within the 
NRQCD factorization formalism. 
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Boyd, Peter Cho, and Ira Rothstein, and the or- 
ganizers of the IV International Conference on 
Hyperon, Charm and Beauty Hadrons. This work 
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ergy under grant number DOE-ER-40682-143. 
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